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Abstract — A framework is developed for analyzing capacity 
gains from user cooperation in slow fading wireless networks 
when the number of nodes (network size) is large. The framework 
is illustrated for the case of a simple multipath-rich Rayleigh 
fading channel model. Both unicasting (one source and one des¬ 
tination) and multicasting (one source and several destinations) 
scenarios are considered. We introduce a meaningful notion 
of Shannon capacity for such systems, evaluate this capacity 
as a function of signal-to-noise ratio (SNR), and develop a 
simple two-phase cooperative network protocol that achieves 
it. We observe that the resulting capacity is the same for 
both unicasting and multicasting, but show that the network 
size required to achieve any target error probability is smaller 
for unicasting than for multicasting. Finally, we introduce the 
notion of a network “scaling exponent” to quantify the rate 
of decay of error probability with network size as a function 
of the targeted fraction of the capacity. This exponent provides 
additional insights to system designers by enabling a liner grain 
comparison of candidate cooperative transmission protocols in 
even moderately sized networks. 

Index Terms — Wireless networking, multicasting, ad-hoc net¬ 
works, sensor networks, cooperative diversity, outage capacity, 
scaling laws. 


I. Introduction 

C OOPERATIVE diversity has been proposed as an at¬ 
tractive approach to combatting slow fading in wireless 
networks [8], [13]. Spatially distributed nodes provide an 
opportunity to create a distributed virtual antenna array and 
can provide substantial gains in slow fading environments. 
There has been a significant interest in studying these gains 
recently; see, e.g., [1], [2], [9], [10], [12] and the references 
therein. 

A convenient channel model for such problems, as has been 
widely adopted in the literature, is a quasistatic one in which 
the parameters are known to receivers, but not to transmitters. 
In such scenarios the classical Shannon capacity is typically 
zero due to the positive probability of the channel experiencing 
an arbitrarily deep fade, so performance is instead quantified 
in terms of outage capacity, which describes the achievable 
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rate subject to a constraint on the level of outage probability 
that can be tolerated [11]. 

Outage analysis applies to a host of multiterminal exten¬ 
sions of such basic channel models as well [14], although 
the expressions become more cumbersome. To address this, 
diversity-multiplexing tradeoff analysis provides a suitably 
coarser scale characterization of such systems by focusing on 
the high signal-to-noise ratio (SNR) regime and examining 
how outage probability scales with SNR in this regime for 
different transmission rates [16]. 

Diversity-multiplexing tradeoff analysis has also proven 
useful analyzing a host of simple network problems. Eor 
example, [15] extends the analysis to the multiple-access 
channel, while [8] extends the analysis to the cooperative 
diversity channel. 

While such analysis of cooperative diversity has proven 
popular, much of the work has been limited to systems in 
which for a given message there is effectively only a single 
destination node and a relatively small number of potential 
relay nodes to participate in the transmission. 

In the present paper, we develop an alternative framework 
within which to examine cooperative protocols. Eirst, our 
emphasis is on the multicasting scenario in which there is 
one message in the network, but generally multiple destination 
nodes. We will focus on two extreme special cases of this 
scenario. One is when all nodes in the network are to receive 
the message, which for convenience we generically refer to 
as multicasting. The other is when exactly one node in the 
network is to receive the message, which we refer to as 
unicasting. 

Second, our framework examines the scenario in which the 
number of nodes in the network is large. This will allow 
us to examine the associated asymptotic scaling behavior of 
cooperative networks. As a by-product, we do not need to 
restrict our attention to high SNR analysis. Indeed, we fix 
the noise power, normalize the channel statistics, and contrain 
the total power transmitted in the network. This allows us to 
parameterize our results in terms of the associated SNR. 

Within this framework, we analyze the relationship between 
transmission rate and the associated error probability. Pro¬ 
vided we use codes of sufficiently long block lengths, outage 
probability dominates the error probability. Specifically, the 
associated outage event is that not all the intended recipients 
are able to decode the message. Not suprisingly, avoiding 
outage in multicasting is more difficult than in unicasting. 

As our main result, we show that under a multipath-rich 
Rayleigh fading network model, a notion of Shannon capacity 
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can be developed. Specifically, there exists a nonzero capacity 
(dependent on SNR) such that for all rates below capacity, 
the error probability can be made arbitrarily small provided 
the network is sufficiently large. Conversely, for all rates 
above capacity, the error probability is bounded away from 
zero regardless of network size. Our achievability result is 
based on a simple two-phase cooperative network protocol 
we develop. By contrast, when one precludes the possibility 
of cooperation, the associated capacity is of course zero. 
Interestingly, our analysis also reveals that despite the fact 
that multicasting outage behavior is dominated by the worst 
node, the multicasting and unicasting capacities thus defined 
are identical. Not surprisingly, we also show that for a fixed 
number of nodes, the probability of error is still much smaller 
in unicasting than in multicasting. 

We further show that finer scale characterizations of be¬ 
havior are possible too. In particular, we define a notion 
of network scaling exponent that characterizes the rate of 
decay of error probability with network size as a function 
of the targeted fraction of capacity. Within this analysis, we 
see, among other insights, that the exponent of our capacity- 
achieving protocol is quite small for rates that exceed half the 
network capacity. 

While our results are specific to our multipath-rich Rayleigh 
fading model, we believe that the associated framework is 
useful more broadly in the analysis of user cooperation gains in 
large networks with more realistic — if more complicated — 
models. Indeed, ultimately our results more generally suggest 
that just as system analysis asymptotic in block length or SNR 
has proven useful, so can one that is asymptotic in network 
size. 

The remainder of the paper is organized as follows. Sec¬ 
tion |II] introduces the network model of interest. The capacity 
result is stated in Section nn It is established by providing a 
converse in Section 113 and an achievability argument in Sec¬ 
tion lYl The scaling of the outage probability with the number 
of nodes is discussed in Section lYll and the network scaling 
exponent is introduced in Section IVTIl Finally, Section IVTIII 
contains some concluding remarks and directions for future 
work. 


II. System Model 

We consider a system with K receiving nodes and one 
source node. For convenience, we label the source node as 
node 0, and the receiving nodes as {1,2,..., K}. Of course, 
in practice, different nodes in the network can act as source 
nodes over orthogonal time or frequency bands as discussed 
in [8]. However, for analysis, it suffices to focus on a single 
configuration. 

We assume a narrowband, slow fading channel passband 
channel model corresponding to a multipath-rich propagation 
environment. In particular, the channel gains hjk between arbi¬ 
trary distinct pairs of nodes (j, k) are independent identically 
distributed (i.i.d.) random variables from a zero-mean, unit- 
variance circularly symmetric complex Gaussian distribution. 
In turn, the signal received at node k at time i is given by 

yk{i)= ^ hjkXj{i),+Zkii), (1) 

ieT(z) 


where T( 2 ) is the set of nodes transmitting at time i, where 
Xj{i) is the symbol transmitted by node j at time i, and where 
Zk{i) denotes circularly symmetric complex i.i.d. Gaussian 
noise of power Nq. Furthermore, the noises among the dif¬ 
ferent receivers are mutually independent. 

In our model, nodes are subject to a half-duplex constraint, 
i.e., a node cannot transmit and receive simultaneously. Thus, 
associated with every valid protocol is a set of binary variables 
of the form Dk{i) that specifies at time i whether node k is 
transmitting {Dk{i) = 1) or receiving {Dkii) = 0). 

The source sends one of M possible messages to the 
destination node(s) over n channel uses (i.e, i = 1,2,... ,n). 
The channel gains between all pairs of nodes remains fixed 
over this duration. In our model, the channel gain hjk is known 
to the receiving node k but not to the transmitting node j. 

We further restrict our attention to protocols in which 
relay nodes cannot revert to receive mode once they begin 
transmitting, i.e., if Dk{i) = 1 for some node k and i < n 
then Dk{j) = 1 for all i < j < n. This restriction precludes 
protocols in which transmitting nodes effectively learn and 
exploit the network channel gains in their encodings. 

Finally, for simplicity, we adopt a long-term sum power 
constraint across the nodes in our model. In particular, with 
Xj{i) for j G 7{i) denoting the (complex-valued) symbol 
being transmitted by node j at time i, we impose an expected 
sum power constraint of the form 


E 


{EE 

i=l fceT(i) 


<p, 


( 2 ) 


where the expectation is taken over the ensemble of channel 
realizations and the set T(i), as well as any other randomized 
aspects of the protocol. Indeed, one will want to consider 
protocols in which the set T( 2 ) depends on the realized channel 
gains. 

The power constraint (13 is a rather natural one for sys¬ 
tems in which there are ergodic channel variations but a 
stringent delay constraint that requires transmission of any 
particular message within a single coherence interval. Indeed, 
although we do not satisfy the sum power constraint during 
the transmission of an individual message, the expectation 
in (|3 ensures that it will be satisfied with high probability 
over a sufficiently long sequence of messages. Nevertheless, 
we remark in advance that the results of the paper do not 
change when the we require the sum power constraint to be 
met with high probability in every coherence interval, and 
our capacity-achieving protocol can be readily extended to 
this case. Ultimately, the expected power constraint merely 
simplifies the exposition. 

The preceding discussion characterizes an admissible pro¬ 
tocol for our analysis, which we formalize in the following 
definition. 

Definition 1: An admissible protocol ttk consists of a set 
of indicator functions {Dk{i)} G {0,1}, which determines 
whether node k is transmitting or receiving at time i; a 
set of encoding functions {fikii)} G C, which determines 
the symbol produced by node k at time r, and a set of 
decoding functions {ipk} G {1;2,...,M}, which determines 
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the message decisions produced by node k at time n. These 
functions are further constrained by their usage as described 
below. 

During the initialization phase of the protocol, the source 
node 0 selects message W G for transmis¬ 

sion. The protocol chooses a priori the sequence Z?o(j) for 
i = l,2,...,n. Without loss of generality, £lo(l) = 1 and 
Dk{l) = 0 for k = 1, 2,..., iT. Moreover, the collections of 
observations at each node k are initialized; = 0. 

At time i, for 1 < i < n, if Z?o(0 = 1’ '^hen source node 
0 uses encoding function 4)o{i) to map W into a transmitted 
complex-valued symbol X[)(i). 

If node k G {1,2,...,iT} is in transmit mode at that time 
(i.e., Dk{i) = 1), the encoding function at node k 

maps y^ and the complex-valued channel gains {hjk,j — 
0,1,..., AT} into a transmitted complex-valued symbol Xk{i), 
which is transmitted over the channel. 

If, instead, the node is in receive mode (i.e., Dk{i) = 0), 
then it collects the complex-valued measurement j/fc(i) and 
updates its set of received symbols via y^ := yfcU{t/fc(z)}. If 
i = n, the decoding function ipk at node k maps y^, and the 
complex-valued channel gains {hjk,j = 0,1,..., iT} into a 
decision Wk- Note that without loss of generality, Dk{n) = 0 
for at least one value of k. If i < n, the node makes a decision 
whether to switch to transmit mode for the remaining duration. 
If it decides to switch, it sets Dk{j) = 1 for i -f 1 < j < n; 
otherwise, it sets Dk{i -f 1) = 0. 

A cooperation-free protocol is a special case of the above 
definition: 

Definition 2: A cooperation-free admissible protocol is one 
for which only the source node transmits, i.e., Dk{i) = 0 for 
1 < k < K and 1 < i < n. 

III. Coding Theorems 

We now develop the relationship between transmission rate 
and error probability for such protocols, in the limit of large 
network sizes. 

We begin with a meaningful definition of capacity. 

Definition 3: A rate R is achievable for the unicasting 
(respectively, multicasting) system if for every network size 
K, there exists an admissible protocol ttk with uk channel 
uses and Mk = 2"^^ messages such that the probability that 
the destination node (respectively, any node) fails to decode 
the message approaches zero as K ^ oo. The (unicasting or 
multicasting) capacity C is the supremum of all achievable 
rates. 

With this definition, we have the following coding theorem, 
which is our main result. 

Theorem 1: The unicasting and multicasting capacities are 
identical and given by' 

C = log(n-A), ,3, 

where P is as defined in (|3 and Nq is the noise power as 
defined via m. 

'Unless otherwise indicated, all logaiithms are base 2. 


Before presenting our proof of this result, we note that to 
achieve capacity — indeed any nonzero rate — requires the 
use of cooperation. Formally, we have the following result. 

Theorem 2: The capacity of cooperation-free admissible 
protocols for both unicasting and multicasting is 

= 0 (4) 

whenever P in 0 is finite and Nq in 0 is nonzero. 

Proof: From Definition|2] a destination node must decode 
directly from the source transmission. Since the channel is 
a Rayleigh fading channel with P/Nq < oo, there exists a 
strictly positive probability of outage — and hence probability 
of error — for every positive transmission rate R. Since these 
probabilities are independent of K, 0 follows. ■ 

We now proceed to the proof of Theorem Q] Since we 
consider large block lengths, error probability is dominated 
by outage probability. Specifically, we consider blocks long 
enough that the probability of error when there is no outage 
is negligible compared to the outage probability. Thus, in the 
error analysis in our proof, we restrict our attention to outage 
probability. 

IV. Prooe of Converse Part 
We develop a converse via a simple upper bound on the 
achievable rate R of Definition |3 In particular, suppose a 
genie conveys the message W to nodes 1,2,...,AT — 1 
and only destination node AT remains to be served. Thus, 
nodes 0,1,..., AT — 1 can coordinate to send the message to 
destination node K. This is clearly a multiple-input single¬ 
output (MISO) antenna system with K antenna elements and 
channel knowledge only at the receiver. Thus, for a given rate 
R, a lower bound on the outage probability for the MISO 
channel is a lower bound on that for both unicasting and 
multicasting systems. We let denote the MISO channel 

outage event. 

To develop such a bound, we first note that it suffices to 
restrict the input distribution to Gaussian. 

Lemma 1 (Teletar [14]): The outage capacity of the slow- 
fading MISO channel with K transmit antennas in 
i.i.d. Rayleigh fading with total power constraint P is 
achieved by an input distribution with a covariance matrix 
diag(Pi, P 2 , ■ • •, Pk), where J2f=i Pj = P- 
We now establish the following lemma; 

Lemma 2: Let e > 0 be arbitrary, and let 

0 ^' 

for a AT-antenna MISO channel. Then the outage probability 
Pr{£^'®®} is bounded away from zero, i.e., 

infPr{£^'®°} > 0. (5) 

Proof: 

Let Pj* ,Pf,... Plf be the power allocations that minimize 
the outage event for the selected rate. The corresponding 
mutual information is given by 


( 6 ) 
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Fig. 1. The two-phase cooperative multicasting protocol. In phase 1 (left), 
the source node (square) broadcasts at a high rate and only a small fraction 
(solid) of the many destination nodes (discs) are able to decode. In phase 
2 (right), these nodes cooperatively broadcast the message to the remaining 
nodes using a suitable space-time code. 


The outage event = {i? > /} is bounded as follows: 


Pr{£^^®°} = 1 


> 1 


> 1 


= 1 


= 1 


Pr(/ > R) 

log 

R 

log 

R 

(} + mT.Lp:) 

log 

R 

^ 

R R' 


In the above derivation 0 follows from the Markov in¬ 
equality, is a consequence of Jensen’s inequality and © 
follows from the Rayleigh model = 1, V£ Since the 

above result holds for all K, 0 follows. 


V. Proof of Forward Part 
A simple two-phase cooperative protocol can achieve any 
rate below the capacity 0, as we now develop. 

A. A Two-Phase Cooperative Protocol 

The protocol of interest is depicted in Fig. ^ Specifically, 
in phase 1, the source node broadcasts the message over ni 
channel uses at a rate i?i and all nodes attempt to decode the 
message. Then, in phase 2, the nodes that are successful in 
decoding the message act as relays and form a virtual antenna 
array, transmitting over the remaining 712 = n — ni channel 
uses at a rate i? 2 - At this point the intended destination(s) 
attempt(s) to decode the message and an outage is declared if 
any of the intended destinations fail. 

Codebook Generation: Suppose that the source generates 
M codewords i.i.d. 62^(0, Pi) for some Pi > 0, each of length 
ni, and all other nodes each generate M codewords i.i.d. 
CAf(0, P2) for some P2 > 0, each of length 77 , 2 . We describe 
the main steps of the protocol for the case of multicasting, 
but indicate the straightforward modifications for the case of 
unicasting. 


Phase 1: The source transmits the codeword correspond¬ 
ing to the intended message from its codebook over ni channel 
uses. We choose the rate in this phase to be (strictly less than, 
but arbitrarily close to^) 

Pi(a) = log2 , (10) 

where 

G(a) = p-i(l-a), (11) 

and P(-) denotes the cumulative distribution function of an 
arbitrary channel gain \hij\'^, and where we have made the 
dependency of Pi on a explicit. 

All nodes attempt to decode the transmission. A node is 
successful in decoding the message if it finds a codeword in 
the source codebook that is jointly typical with the received 
sequence. Let Ki be the number of nodes that are successful 
in decoding the message from the source. We label these nodes 
as 1,2,..., Ki. These nodes participate in phase 2. 

Phase 2: Each of the Ki nodes successful in decoding 
the phase 1 transmission next transmits the corresponding 
codeword from its codebook over 772 channel uses. The rate 
in this phase is set to (strictly less than, but arbitrarily close 
to) 

P2(/3) = log (^1 + ^(1-/?)) , (12) 

where 0 < /3 < 1 is design parameter, the dependence of P 2 
on which we have made explicit. 

Each of the remaining K 2 = K — Ki nodes attempts 
to decode the message at the end of the second phase. 
Node k, upon receiving is observations^ finds a message 
Wk and a subset of nodes Sfc = {A:i, fc 2 ,..., A:|§j.|} C 
{1, 2,..., K} \ {fc} such that the corresponding set of code¬ 
words {x'^^{wk),x'^^{wk), ■. ■ is jointly typical 

with 7 /^^. It declares the message Wk to be the transmitted 
message if a unique pair {wk, §k) exists and declares a failure 
otherwise. 

In the case of unicasting, if the destination node is suc¬ 
cessful in decoding the message in phase 1 then it does not 
participate in phase 2. Otherwise it continues to listen to the 
transmissions and attempts to decode the phase 2 transmission. 
An error occurs if the destination fails to decode the phase 2 
transmission. 


B. Protocol Analysis 

Eirst, we analyze code rate. To begin, it is straightforward to 
verify that 771 and 772 are completely determined by the choice 
of rates. In particular, let the overall rate of our protocol be 
R, so that there are M = 2”^ possible messages to send over 
77 channel uses in the system. Then it follows that 


mRi = 772 i ?2 = logM. (13) 

^This technicality ensures that the probability of error when not in outage 
will approach zero uniformly over all channel realizations. A similar techni¬ 
cality applies to the rate in phase 2. 

^We use the superscript ^^2 to denote the vector formed from the 712 
variables corresponding to time instants ni + 1, ni + 2,..., n, i.e., phase 2 
of the protocol. 
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From o it is straightforward to calculate the overall 
effective rate R of the system. In particular, since log M = nR 
and since ni + n 2 = n, J13> implies that R satisfies 

R ^ Ri{a) ^ i?2(/3)' 

Second, we analyze the power constraint 0 . In phase 1, 
the transmitted power is Pi, so provided Pi < P, our power 
constraint is met in this phase. To analyze the power used 
in phase 2, we begin by noting that on average a fraction a 
of the nodes are able to decode the message after phase 1. 
Specifically, the number of nodes Ki successful in phase 1 
has mean 

E[Ki] = aK (15) 

since Ki is a binomial random variable, viz., (cf. (HOj) 

K 

^1 = ^ l{|/iOiP>G(a)}, (16) 

i=l 

where Ipj is an indicator function, which equals 1 if its 
subscript is true and 0 otherwise, and where we have set 
Pi = P. Hence, 


The following lemma provides bound on the conditional 
outage probability that will be useful in the sequel. 

Lemma 3: Suppose ki > aK{l—(3) nodes are successful at 
the end of phase 1. Then the conditional probability of outage 
is given by 


Pr{£ 


UC 

a,l3,K 


Ki = ki} < 




(19) 


Proof: To obtain d, it suffices to bound 

Py{Gk{Ki) < I - P \ Ki = ki} since Pr{£“^_^ | 
Ki = ki\ = Vr{GK{Ki) < 1 — /3 | Ki = fci}. This can be 
accomplished by the Chernoff bound, since in (II 8> 

can be written as 


fci 

GK{ki)=Y.g, 


i=i 


where the 



( 20 ) 


( 21 ) 


E 




j=i 


= P2E[Ki] = P2aK, 


from which we see that the power constraint is satisfied in 
phase 2 provided P 2 < P/aK. 

Finally, we analyze the outage probility, i.e., the probability 
of outage of a node that is unable to decode at the end of the 
protocol. For convenience, let us exploit symmetry and label 
this node K, while the nodes that are successful in phase 1 we 
label 1, 2,... iTi. From straightforward MISO system analysis, 
node K will fail to decode the message whenever Ki < K 
and 

i?2(/3)>log(^l + GK(iTi)^j, (17) 

where ^ 

= (18) 

f=i 

is the effective MISO channel gain of node k, and where we 
have set P 2 = P/{aK). But since R 2 was chosen according 
to lll2> in phase 2, (I17> implies that outage will occur when 
Gk{Ki) <1-/3. Accordingly, the outage events £™^ ^ and 
£™^ if for unicasting and multicasting, respectively, take the 
form 


if-i 

cue _ I I /lUC 

^(y,f3,K '^ki 

ki—0 
K-l 

cmc _ I I /I m.( 

^oc,f3,K — '^ki 

ki—0 


where, for ki £ {1,2,..., K}, 

All = {Ki = ki, Gif(/ci) < 1-/3} 
AT:={Ki = ki, ^^mm^G.(fc 0 <l -/3 


are i.i.d. random variables. Specifically, we obtain, for any 
s > 0, 


Pr{Gif(iTi) < 1 -/3 I ifi = A:i} 


= Pr 


l^g-sGK(ifi) 







g-sGK(ifl) 



gs(l-/3) 

“ (1 + ’ 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


where (EUi follows from the Markov inequality, (d follows 
by from ( I20> . and ( I26l l follows from evaluating the character¬ 
istic function of the exponential random variables (ED- 
In turn, since (ED holds for all s > 0, we can choose the 
particular value 

Substituting \21\ into (I26> yields (I19> as desired. ■ 

We now show that the probability of outage can be made 
arbitrarily small by selecting aK appropriately. 

Proposition 1: The probability of outage in unicasting de¬ 
creases exponentially with aK for every /3 > 0. Specifically, 
for every 0 < e < /3, 


Pr{£“/3.if} = E [Pr{£“^,^|iFi}] 

< exp(—aiTe^/4) -f exp {aK^{f3, e)} , 


(28) 


where 


7(/3,e) =/3-e-f (1 



< 0. (29) 













6 


Proof: 

(0,/3), 


To obtain (HHJl we observe that, for some e G 


Pr{£ 


UC 

a,l3,K 


} 


E I = fci} Pr{iTi = fci} 

ki:ki<aK{l — e) 

+ E I K, = fci} Pr{iTi = k,} 

ki:ki>aK(l — e) 


(30) 

< PrjiTi < Q!itr(l — e)} + max Pr{£“/ 3 ,/s- | Ki = ki} 

ki>aK{l-e) 

= PrjiTi < aK{l - e)} + Pr{£“^ k \ Ki = aK{l - e)}, 

(31) 


where (13 U exploits that outage probability is a decreasing 
function of ki. Finally, using the binomial Chernoff bound 
(see, e.g., [4]) 

PrjiTi < aK{l - e)} < 


for the first term in (131> . and applying Lemmal^to the second 
term, yields (ED as desired. ■ 

In turn. Proposition [2 can be used to bound the correspond¬ 
ing probability of outage in multicasting. 

Proposition 2: The probability of outage in multicasting 
decays asymptotically with aK. Specifically, 


^ [P>-{£r/3.ifl^i}] < if Pr{£“^,^}. (32) 
Proof: First, we bound the conditional outage probability 
according to 

Pr{£“^^_^|ifi = fci} (33) 


= Pr{ 

min 

GkiKG <l-3\Ki = kA 

•c:ki ‘^k'^K. 



f ^ 



= Pr ■ 

u 

{G^iKi) <1-/3} 

II 


i=Ki + l 


<{K-ki) Pr{£“^_^|ifi = fci} (34) 

<ifPr{£^=^_^|ifi = fci}. (35) 

where (ED is a simple application of the union bound. Taking 
the expectation of both sides of (E3 with respect to Ki, 
we obtain (l.32t . Finally, since Proposition [Q establishes that 
the unicasting outage probability decays exponentially, (E3 
implies that the multicasting outage probability does as well. 


Propositions E] and El can be used to establish the forward 
part of the coding theorem for both unicasting and multicast¬ 
ing. 

Proof of Theorem El To show that our two-phase protocol 
can approach the capacity we show that the outage proba¬ 
bility can be made arbitrarily small while operating arbitrarily 
close to the capacity. Suppose that a > 0 and /3 > 0 are 
arbitrary. Since the outage probability decreases exponentially 
in K, we can choose a K large enough to make the outage 
probability sufficiently small. Next, note that by choosing a 
and (3 sufficiently small, we can make Ri{a) sufficiently large 
and i?2(/3) sufficiently close to C. As particular examples, it 
suffices to take a ~ 1/ log if and (3 ~ 1/if so that a,P ^ 0 


but aK oo. Thus, we can have the effective rate (fT?l be 
arbitrarily close to C, while keeping the outage probability 
sufficiently small. 

■ 

An intuition behind the achievability result is that in the 
limit of a large number of nodes, we can find sufficiently many 
nodes (albeit a small fraction of the population) with very large 
channel gains and they can be served over a small number 
of channel uses in the first phase (i.e., ni is a negligible 
fraction of 712 ). These nodes then simultaneously cooperate 
to serve the remaining nodes. Since sufficiently many nodes 
are transmitting in the second phase, we have enough diversity 
in the system to drive the outage probability to zero. 


C. Multiple Antenna Generalization 

It is possible to generalize our results to the case where 
the at least some of the nodes in the network have multiple 
antennas. In particular suppose that the node i has Ti antennas. 
In the case of unicasting, our two-phase protocol can be 
straightforwardly extended to obtain the following; 

C = T,,log(l + A). (36) 

In the bound (El the key quantity of interest is Tk, the 
number of antennas at the destination node; the number of 
antennas at the source and relay nodes do not impact capacity. 
Note that, in the first phase of the protocol, we can still 
communicate to a large number of relay nodes, regardless 
of the number of antennas at each relay. These nodes then 
form a virtual antenna to communicate to the destination in 
phase 2. This reduces to the case of a multiple-input multiple- 
output (MIMO) system when the number of transmit antennas 
is much larger than the receive antennas. Using the channel 
hardening result for such systems — see, e.g., [7] — one 
can establish that rate is achievable. The converse is 

analogous to the single antenna case in Section CYl 

An analogous argument for multicasting can also be devel¬ 
oped, from which we have the following: 

C = min{Ti, T 2 , ...,Tk}- log (l + ^) • (37) 

Evidently, (ED can be much smaller than (ED — the lower 
bound system rate is limited governed by the node with the 
fewest antennas in multicasting rather than the destination 
node. 


VI. Outage Scaling Behavior 

Our capacity result determines the rates for which outage 
probability goes to zero with increasing network size for multi¬ 
casting and unicasting. Often, a finer grain analysis is required 
by system designers. In this section, we develop the manner in 
which outage probability goes to zero with increasing network 
size for the two-phase protocol of Section lYl which provides 
several additional insights. For example, while we have shown 
that multicasting and unicasting share the same capacity, here 
we show how their respective outage probability curves differ. 
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A. Outage Probability Approximations 

While J28> and (I32> bound the outage probabilities of 
interest, these bounds are not tight. Nevertheless, good ap¬ 
proximations to the actual outage are readily obtained, as we 
now develop. 

The outage probability of a unicasting system under the 
two-phase protocol can be approximated by 


~ ■^exp{-aiT(l - /?)} X 


where 


with 


exp < —KD{'y\\a) + yKln 


7 = 


t(l -/3)e 


P = 



(38) 

VI+ 4/7-1 

(39) 

2 /i 

aHl-P) 

1 — a 

(40) 


In turn, the outage probability of the multicasting system 
can be approximated in terms of this unicasting approximation 
according to 


= 1 - (1 - « if Pr{£“^,;,}. 

(41) 

A derivation of the approximation (El is provided in the 
Appendix. 


B. Accuracy of Outage Probability Approximations 

In this section, we compare our outage probability bounds 
(I28> and (I32> : and our approximations (138> and (141> . to the 
actual probabilities via Monte Carlo simulations. In particular, 
we choose a target rate below capacity and evaluate the outage 
probability as a function of the network size K. We evaluate 
the expectations over Ki in the bounds (El and E3 by 
numerical integration. 

For our comparison, we set a rate of 

i? = i log f 1 -f —) , 

2 ^ V 2No) 

which is 1/2 of capacity in the high SNR regime and 1/4 of 
capacity in the low SNR regime. This rate point is realized 
by the parameter settings G{a) = 1/2 and /3 = 1/2 in our 
two-phase protocol, so i?i(a) = i?2(/3)) in (I10> and (ED, 
respectively. 

Fig-El depicts the results. Several observations are worth 
emphasizing. 

Remarks: 

1) First, the outage curves for both unicasting and mul¬ 
ticasting approach zero with our cooperative protocol, 
which is a consequence of the transmission rate being 
below capacity. Note that, by contrast, for cooperation- 
free admissible protocols, the outage curves will not 
decay with network size. 

2) Multicasting incurs significant penalty over unicasting 
in terms of outage probability for a fixed network size 
K. In particular. Fig. El confirms that the multicasting 



Fig. 2. Outage probability for unicasting and multicasting via our two- 
phase protocol as a function of network size. The solid curves correspond 
to our Monte Carlo simulations, the dashed-dotted curves to our analytical 
approximations, and the dashed curves to our bounds. The top set of curves is 
for multicasting; the bottom set for unicasting. In the protocol we set (3 = 0.5 
and i?i(a) = R, so that the rate is R = (1/2) log(l + P/{2Nq)) < C. 

outage probability is indeed roughly a factor K larger 
than the unicasting outage probability. 

3) The slopes of the outage log-probability curves are 
asymptotically constant, and the bounds are good pre¬ 
dictors of the asymptotic slopes. This is perhaps not 
surprising since we used Chernoff techniques to derive 
the bounds. Indeed, in many communication problems 
the Chernoff exponent is close to the correct exponent. 
However, the bounds are not particularly close to the the 
correct outage curves. 

4) The analytical outage probability approximations are 
asymptotically quite close to the true curves, converging 
to within a factor of roughly 3 in probability for large 
network sizes. In addition, these approximations appear 
to be actual upper bounds at least in case study depicted, 
though this is conjecture. 

5) The asymptotic slopes of the outage log-probability 

curves for both unicasting and multicasting are identical. 
In the next section, we will develop this slope as the 
network scaling exponent of the protocol, which we 
denote using For a target outage level, this slope 
can be used to quantify the asymptotic network size 
gap between unicasting and multicasting. In particular, 
suppose that for a fixed choice of a and (3 in the pro¬ 
tocol, nodes are required to achieve some target 

outage probability e in unicasting. Then the number of 
nodes required to achieve the same outage probability 
in multicasting is, asymptotically, 

(e) = ^ log (42) 

Ens 

To verify (l4^ . it suffices to recognize that the vertical 
distance between the unicasting and multicasting outage 
probabilities is, in accordance with (EB, asymptotically, 
logA:“'^(e). 







VII. Network Scaling Exponent 


In this section, we explore, in more detail, the asymp¬ 
totic rate of decay of the outage probability with network 
size, which we have termed the network scaling exponent. 
This exponent captures meaningful information for system 
designers. For example, at the transmission rates to which 
Fig. 121 corresponds, outage probabilities for our two-phase 
protocol decay reasonably quickly in a practical sense — i.e., 
the network scaling exponent is reasonably large. However, 
as we will see, at rates close to capacity, it turns out that 
outage probabilities decay very slowly as a function of network 
size, corresponding to a small network scaling exponent. This 
implies that very large network sizes may be needed to achieve 
practical target error rates. 

Before beginning our development, note that the network 
scaling exponent is the natural counterpart to the classical 
error exponent for traditional channel codes. In particular, the 
classical error exponent captures the exponential rate of decay 
of error probability with block length as a function of the 
the targeted fraction of capacity; see, e.g., [6]. Analogously, 
the network error exponent captures the exponential rate of 
decay of error probability in unicasting and multicasting with 
network size as a function of the targeted fraction of capacity. 

Formal definitions follow. 

Definition 4: The network reliability function with respect 
to a sequence of admissible protocols t^k in Dehnition ^ is 
given by 

£,;({«}) = - l.m ,43) 

where denotes the outage event for a system with K 
nodes under the protocol -kk- 

Definition 5: The network scaling exponent is the supre- 
mum of the network reliability functions of all sequences of 
admissible protocols with a rate that is at least a fraction r of 
the capacity at a given SNR, i.e., 

F;ns(r, SNR)= sup S” ({7r_R-}), (44) 

{iTK-I sTR.SNR) 

where J’(r, SNR) is a set of sequences of admissible protocols 
with a rate that is a fraction r of the capacity. 

The following establishes that, as with capacity, unicasting 
and multicasting are not distinguished by their network scaling 
exponents. 

Proposition 3: The network scaling exponent is the same 
for both unicasting and multicasting. 

Proof: First, for any sequence of admissible protocols, 

Pr{£“^}>Pr{^-}, 

so that 

EZ{r, SNR) < SNR). (45) 

Furthermore, if {tt^} achieves the supremum for the uni¬ 
casting system, then from a simple application of the union 
bound it follows that, for each K, 

Pr{£“f}<ifPr{£-}, 



0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 


fraction r of capacity 

Fig. 3. The upper envelope of the plotted points indicates the network scaling 
exponent for the two-phase cooperative protocol as a function of the targetted 
fraction of capacity r. Each point corresponds to a particular value of a and 
/3. In this example, SNR = 0 dB. 


and hence 


E^%r, SNR) > Efrm^K}) = KsiP SNR). (46) 


Combining §5} and (l4^ we obtain SNR) = 

SNR) as desired. ■ 

In the remainder of this section, we analyze a lower bound 
on the network scaling exponent by optimizing over the class 
of the two-phase protocols described in Section 0 For a hxed 
choice of a and we can express the fraction of the capacity 
achieved by the protocol as 


r{a, /3, SNR) 


R{a,/3, SNR) 
C(SNR) 


(47) 


where we have made the dependency of both R and C in (HI 
and 0 . on the parameters of interest explicit. We define the 
network reliability function of the user cooperation protocol 
in Section m as 


Fins (r, SNR) = sup 

a,/3:r(a,/3,SNR)<r 


lnPr{£a,/3,iT} 

R —>00 K 


(48) 


which constitutes a lower bound on Ens{r, SNR) in il44t . Note 
that in the above definition, we have constrained a and (3 to 
be constants independent of K. 

The upper envelope of the points in Fig.|3]indicates a lower 
bound on the network reliability function of our two-phase 
protocol. Each point in the plot corresponds to a particular 
choice of a and /3 in the protocol, for which we have 
numerically evaluated Ef^ in (Hill for different values of r 
at SNR = 0 dB. 

Perhaps the most striking observation from Fig.|3lis that the 
error exponent for the two-phase protocol is quite small when 
aiming for rates that are more than about half of capacity. 
This implies that while the protocol is capacity achieving, 
it may require a prohibitively large number of nodes to 
achieve rates anywhere close to this capacity. It remains to be 
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determined whether there exist more sophisticated protocols 
with substantially higher exponents in this regime. 

As a final comment, it should also be noted that Fig. 0 
effectively characterizes the efficient operating frontier for the 
protocol. In particular, given a network with K nodes and an 
allowable outage probability e, one can approximate the scal¬ 
ing exponent by —hiejK and determine the corresponding 
value of r, which is an estimate of how close one can expect 
to get to capacity in the system. 


VIII. Concluding Remarks 


Perhaps the main contribution of this paper is a framework 
for analyzing user cooperation protocols in the limit of large 
network size (number of nodes), which we have illustrated 
in the case of a multipath-rich Rayleigh fading environment. 
Within this framework, we have introduced a meaningful 
notion of Shannon capacity for this regime and presented a 
simple two-phase protocol that can achieve rates arbitrarily 
close to capacity. A finer grain analysis of this two-phase 
protocol in terms of its network scaling exponent, which 
characterizes the rate of decay of error probabilty with network 
size, shows that it may require prohibitively large number of 
nodes to achieve rates close to the capacity with this protocol. 

One important direction of future work is to study more 
sophisticated models beyond the Rayleigh fading model within 
our framework. One could for example incorporate the effects 
of network geometry and shadowing into the model. More 
generally, it would be of interest to study a class of channel 
models for which user cooperation plays a fundamental role in 
enabling reliable communication in multicasting. The Rayleigh 
fading model considered here clearly belongs to this class, but 
we believe the class may be quite rich and may include many 
other models of practical importance. 

Another important direction is to investigate how system 
performance changes when the sum power constraint is re¬ 
placed with individual power constraints. With individual 
power constraints, the system capacity will increase with the 
number of nodes — in fact, the MISO upper bound increases 
according to 0(logiT). It remains to be determined whether 
there exist cooperative multicasting protocols that approach 
this upper bound or whether one can develop tighter upper 
bounds for this scenario. 

Finally, as noted in Section IWI the two-phase protocol 
may require prohibitively large number of nodes to achieve 
rates close to the capacity. It remains to investigate whether 
more sophisticated protocols can improve the network scaling 
exponent substantially in this regime. 


Appendix 

Derivation of Outage Approximation 

First, we write Pr{£“^ jf} form 

K 

= ^ Pr{ifi = fci} Pr{£-^,^ I K, = k,} 

k-y — l 

K 

« ^ Pr{iTi = k,} Pr{£-^,^ | = k,}. 

ki=a(l-l3)K 

(49) 


Note that we have dropped the contribution of terms with ki < 
a(l — (3)K in the summation, since we expect their aggregate 
sum to be small as they deviate significantly from the mean 
E[Ki] = aK. 

We now approximate each of the two factors in J49t . The 
right factor we approximate by the upper bound ( I19t . The 
left factor we replace with via Stirling’s approximation for 
binomial distributions [5, p. 284], yielding 


Pt{Ki 


fci}= 


1 

-^= exp 




(50) 


where T>(-|j-) denotes the binary relative entropy function, i.e., 
for any 0 < p,q < 1, 


D{p\\q) =p\n^ + (1 -P)ln 

and where a is the parameter of Ki (cf. ( I16H . 
Thus, substituting J19t and (EOj into J49t yields 

Pr{£:i!/3.A} 




ki=a(l-/3)K 


aK{l — j3)e 
ki 


(51) 


K 


exp {—aAr(l —/?)}. (52) 


Finally, we approximate by an approximation to the 
largest single term in the summation, viz.. 


Pi'{£“/3.a} ~ ■^exp{-aA:(l - /))} x 

max eyipi—KD{'y\\a)-\- 
7e(a(l-/3).l) i V 7 

(53) 

Since the term in the exponent being minimized in o is 
differentiable and convex in 7, the optimizing 7 is the value 
at which the associated derivative is zero, i.e., 

7 ^ 

= F (54) 

1-7 

where p is as given in ( I40> . Finally it is straightforward to 
verify that J54> has a solution in (a(l — (3), 1) and it may be 
solved explicitly, yielding ( I39> . 
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